function [out] = geometric_mean(x)

x = x(:);

x = nonzeros(x);

log_x = log(x);

out = exp(mean(log_x));

end

